Phys. Rev. B, in press 



cond-mat/0406431 

Effects of the On-Site Coulomb Repulsion in Double Exchange Magnets 
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We investigate the zero-temperature phase diagram and spin-wave properties of a double exchange magnet 
with on-site Hubbard repulsion. It is shown that even within a simple Hartree - Fock approach this inter- 
action (which is often omitted in theoretical treatments) leads to qualitatively important effects which are 
highly relevant in the context of experimental data for the colossal magnetoresistance compounds. These 
include the asymmetry of the doping dependence of spin stiffness, and the zone-boundary "softening" of 
spin wave dispersion. Effects of Hubbard repulsion on phase separation are analyzed as well. We also show 
that in the ferromagnetic phase, an unusual temperature-dependent effective electron-electron interaction 
arises at finite T. The mean-field scheme, however, does not yield the experimentally observed density of 
states depletion near the Fermi level. We speculate that proper treatment of electron-electron interactions 
may be necessary for understanding both this important feature and more generally the physics of colossal 
magnetoresistance phenomenon. 

PACS numbers: 75.47.Gk, 75.47.Lx, 75.30.Ds, 75.10.Lp 



I. INTRODUCTION 

The phenomenon of colossal magnetoresistance (CMR) 
is known to occur in a broad group of compounds, corre- 
sponding to different crystal structures, chemical compo- 
sitions, and doping levels [1]. In addition to various heav- 
ily doped manganese oxides, the CMR effect is also ob- 
served in certain magnetic semiconductors and spinels [2] ; 
it is natural to expect that in all these cases, the physical 
origins of the CMR are similar. Thus, a proper minimal 
theoretical model of a CMR system should account for 
the important common features shared by all these ma- 
terials, while leaving out the peculiarities of crystal envi- 
ronment and atomic structure of individual compounds. 
It is universally recognised that one of these common fea- 
tures is the sizable ferromagnetic Hund's rule coupling, 
Jh , between the spins of magnetic ions and those of con- 
duction electrons, which gives rise to the double exchange 
ferromagnetism [3] of the CMR compounds. The purpose 
of the present article is to draw attention to the fact that 
another ubiquitous intra-atomic interaction, namely the 
Coulomb (Hubbard) repulsion U, also affects magnetic, 
electronic, and transport properties of the system in a 
profound way, and may play a crucial role in the basic 
physics of the CMR. While some effects of this interaction 
have been addressed in the past (see, e.g., Refs. [4-7]), its 
potential importance is not yet fully appreciated. We will 
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argue that the on-site Coulomb repulsion strongly affects 
the magnetic properties of the system; some generic ex- 
perimental facts are recovered. We also suggest that the 
effects of Hubbard repulsion merit further investigation 
beyond the mean-field approach. 

We start with the standard double exchange Hamilto- 
nian, supplemented with a Hubbard repulsion term: 
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Here, the fermionic operators Cja correspond to con- 
duction electrons, hopping between the atomic sites of 
magnetic ions with spins Si, and the vector a"'' is com- 
posed of Pauli matrices. Electron concentration is de- 
noted by X, hence the hole density is given by 1 — 2: 
(we note that in the experimental literature on the CMR 
manganates, the opposite convention is often used). In 
order to discuss our results within the context of experi- 
mentally observed magnetic phase diagrams of the CMR 
manganates, we also include the antiferromagnetic su- 
perexchange coupling J between the ionic spins. We 
will treat the ionic spins as classical, S* ^ 1; quan- 
tum corrections are not expected to modify the effects 
of Hubbard term in a qualitative way. For the case of 
CMR manganates, S = 3/2, and the band theory calcu- 
lations [5,8] suggest the typical values of i '--^ 0.3 — 0.5 
eV, Jh ^ 2.5 eV, and U ^ 8 eV. The value of J can 
be roughly estimated from the experimentally observed 
Neel temperatures in the fully doped (no conduction eg 
electrons, x = 0) case [9-11], T/v -- 100- 200K, yield- 
ing J ^ 5 — lOmeV. We will consider the case of square 



(2D) or simple cubic (3D) lattice, assuming that the lat- 
tice spacing is equal to unity. Throughout the paper, 
chemical potential is denoted by ;U — (J/f /2). 

While the important and highly non-trivial effects of 
the orbital degree of freedom in the CMR manganates 
are of great interest to both theorists and experimental- 
ists working in the field [1], in writing Eq. (1) we assumed 
that there is only one atomic orbital available to conduc- 
tion electrons at each site. The reasons for this drastic 
simplification are three-fold: (i) in the CMR manganates, 
orbital structure is strongly dependent on the crystalline 
environment and varies for different compounds and dop- 
ing levels [12]; (ii) yet another situation takes place for 
magnetic semiconductors exhibiting CMR, like EuS or 
EuSe, where the three t2g conduction bands show no 
Jahn-Teller splitting; therefore at present it seems dif- 
ficult to conclude that a realistic treatment of orbital 
effects is crucial for understanding the basic physics of 
the CMR phenomenon; and (iii) we will see that the ef- 
fects of Hubbard repulsion are rather complex already 
in the single-orbital case; we believe that these should 
be understood before a more complicated model [6] is 
advanced. 

We begin with a brief overview of the low-temperature 
properties of the double exchange model, Eq. (1) with 
U = 0. The carrier spectrum in the ferromagnetic state 
of a 2D (3D) system is given by el'^ = =f(Jh/2) + e^: 
with ej: = — tcosfca, — t cosky{—t coskz) ■ A sufficiently 
large value of Hund's rule coupling Jh then results in 
a complete spin polarisation (half-metal) of conduction 
electrons within the entire doping range of < a; < 1, 
in agreement with experimental data [13] (see, however, 
Ref. [14]). At Jh -^ oo and J = 0, any deviation of 
ionic spins from ferromagnetic order results, via the dou- 
ble exchange mechanism [3] , in a narrowing of the spin- 
up conduction band and therefore costs positive energy. 
The corresponding value of spin stiffness D is then pro- 
portional to conduction band energy, DS — \E\/4:d with 
E ~ J epri:j^df'k/{2'!TY (where n^: is the Fermi distribution 
function and d is the dimensionality of the system) and 
is therefore symmetric in electron density x with respect 
to the quarter-filling, x = 1/2. In a more general case 
of finite Jh and J > 0, this double exchange ferromag- 
netism competes against antiferromagnetic tendencies, 
which originate from two distinct physical sources. In ad- 
dition to the direct superexchange contribution J (which 
is responsible for the antiferromagnetism of the system 
at X = and, roughly, can be assumed to be doping- 
independcnt), there arises an indirect antiferromagnetic 
interaction [15] which further lowers the relative energy of 
antiferromagnetic phases. This interaction is due to vir- 
tual transitions of conduction electrons between the two 
components of the spin-split band [16], and its strength 
increases with increasing electron density x. Indeed, the 
net antiferromagnetic contribution to the spin stiffness 



DS of a double exchange ferromagnet at x ^ 1 equals 
-J-TTxHy{2JH) in 2D and - J- (67r)2/3x5/3t2/(10Jff) 
in 3D, and grows to —J — t^ /{4,Jh) at a; = 1. While 
the actual destabilisation of the ferromagnetic phase with 
increasing strength of antiferromagnetism proceeds via 
phase separation [2,17,18], rather than a spin-wave me- 
diated phase transition, this behaviour of spin stiffness is 
in line with the overall conclusion on the phase diagram 
asymmetry: in the U = case, antiferromagnetic tenden- 
cies are more pronounced at x > 1/2 than at x < 1/2. 

This expected behaviour does not agree [4] with the ex- 
perimentally observed low-temperature magnetic proper- 
ties of the CMR manganates [9-11]. In broad terms, it 
is fair to say that the CMR manganates are ferromag- 
netic [19] at x > 1/2 and antiferromagnetic at x < 1/2. 
Although the presence of narrow ferromagnetic regions 
(or possibly ferro-antiferromagnetic phase separation) at 
X < 1/2 has been reported in some cases, the wide fer- 
romagnetic area is always located at low hole doping, 
X > 1/2. The investigation of lightly-doped manganates 
with 1 — a; ^ 0.1 is complicated by the sample prepara- 
tion issues. So far, only the 3D perovskite materials are 
available in this region; these are typically [9,10] found 
to remain ferromagnetic down to the very low values of 
1 — X, with a likely exception of the x = 1 endpoint [20]. 
This is in contrast with robust Neel antiferromagnetic 
ordering, characteristic for all manganates at a; <C 1. 

This qualitative discrepancy can be alleviated by tak- 
ing into account the on-site Coulomb interaction, U. The 
latter does not affect the energy nor the carrier spec- 
trum of the fully-magnetised half-metallic ferromagnetic 
state; however, when the neighbouring spins (say. Si and 
Sj) are out of alignment, there arises a non-zero hop- 
ping matrix element [21], connecting the conduction elec- 
tron state at site i with the electron spin directed along 
Si, and the state at site j with the spin anti-aligned 
with Sj. Thus, when two electrons are placed on sites 
i and j, there is a non-zero quantum-mechanical prob- 
ability of double occupancy on-site, and the associated 
Coulomb energy: the electrons repel each other. There- 
fore, Hubbard interaction enhances ferromagnetism, and 
the strength of this effect increases with x [22]. More- 
over, at a finite temperature T, when the ionic spins are 
misaligned due to thermal fluctuations, an unusual T- 
dependent electron- electron repulsion arises in the ferro- 
magnetic phase (in Sect. Ill below we will see that there 
also arises another, essentially many-body, contribution 
to the effective electron-electron interaction). We note 
that both of these effects, which we will consider in some 
detail below, are absent in the widespread simplified pic- 
ture of double exchange, when the value of Hund's rule 
coupling is assumed to be infinite, making the double 
occupancy impossible. 

Although U is, in fact, the largest energy scale in the 
problem, we will use the Hartree-Fock approximation. 



which formally becomes accurate only at small values of 
U/ Jh- It is, however, expected that substituting larger 
values of U into our equations should yield the estimates 
which are adequate at the qualitative level. We note that 
in general, this "stretching" of the Hartree-Fock scheme 
requires some caution, as, for example, the energy of a 
single spin-down particle in the ferromagnetic state and 
in the presence of a partially filled spin- up band of width 
W -^U IS clearly of the order of J^ + T4^ [23] . This is in 
contrast to the Hartree-Fock result for the energy of the 
spin down electron. 



ei-4 + '^'5' Js^Jh+xU 
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[here, Js is the mean- field (Stoner) band splitting]. How- 
ever, the contribution of the Coulomb energy to the prop- 
erties studied in the present paper originates from an 
integral over many intcr-subband contributions (or over 
many spin-down electron states); it is hoped that the 
Hartree-Fock mean field-type approximation is more re- 
liable in such a case. 

We will be interested in the experimentally relevant 
case when the value of Js is large in comparison to the 
Fermi energy, ep- It should be noted that in this regime, 
the family of models with the Hamiltonian (1) and differ- 
ent values of the ratio xU / Jh provides a connexion be- 
tween the conventional double exchange system [U = 0) 
and the large-C/ Hubbard model {Jh ^0). In fact, due 
to the considerable uncertainty in the values of t and 
Jh quoted in the literature, it is not clear whether Jh 
alone could always account for a complete carrier spin- 
polarisation in the ferromagnetic state of the CMR man- 
ganates [14]. It is thus possible that in real systems, the 
half-metallic state (which at the mean-field level is im- 
plied by the condition Js > ef) would not have been 
reached without further enhancement of band splitting 
by the on-site Coulomb repulsion U [cf. Eq. (2)]. If 
this is indeed the case, it might lead to potentially im- 
portant and novel many-body effects, at both zero and 
finite temperatures. These lie beyond the mean-field ap- 
proach taken in the present work - within the Hartree- 
Fock scheme based on Eq. (2), it is indeed unimportant 
whether the perceived half-mctallicity is partly due to 
the effects of Hubbard repulsion. 

The spin wave theory of double exchange ferromag- 
nets in the presence of on-site Coulomb repulsion is con- 
structed in Section II. We evaluate the spin stiffness, 
D, and show that the on-site interaction strengthens fer- 
romagnetism, while restoring the correct asymmetry in 
the doping dependence of D. In addition, we show that 
this interaction results in a suppression of magnon energy 
near the zone boundary^ in comparison with the nearest- 
neighbour Heisenberg dispersion law. We then calcu- 
late the strength of the novel temperature-dependent 
electron-electron interaction (Sect. III). While this in- 
teraction appears negligible in the manganates, in case of 



lightly doped CMR magnetic semiconductors it does lead 
to an appreciable renormalisation of nearest-neighbour 
Coulomb repulsion and (correlated) carrier hopping am- 
plitude. 

The effects of Hubbard repulsion on phase separation 
in double exchange magnets at T = are discussed in 
Sect. IV. In the C/ = case, the zero-temperature 
phase diagram (much like the doping dependence of spin- 
stiffness) suggests that the area of stability of the ho- 
mogeneous ferromagnetic state is shifted towards the 
electron-doped end, x < 0.5, which is at variance with 
generic experimental observations (see above). We show 
that inclusion of U alleviates this difficulty as well. 

The observed suppression of the carrier density of 
states near the Fermi level at low temperatures [24] is 
likely to be of the same origin as the much more pro- 
nounced depletion of the density of states [25,26] (some- 
times termed "pseudogap") in the vicinity of the Curie 
temperature. It appears, in turn, that in order to ad- 
equately describe the physics of CMR phenomenon one 
has to understand the nature of the pseudogap. It is 
therefore important that a proper description of a low- 
temperature ferromagnetic state of the CMR manganates 
should include the correct energy dependence of the 
density of states. In Sect. V we show that standard 
Altshuler-Aronov mechanism utilizing a combination of 
Coulomb repulsion (or the effective electron-electron in- 
teraction derived in Sect. HI) and impurity scattering 
cannot account for the measured depletion of the density 
of states. This signals the insufficiency of our mean-field 
treatment in this case, and the presence of strong energy- 
dependent correlation effects even at low temperatures. 

The implications of our findings are further sum- 
marised in Sect. VI, where we also discuss prospective di- 
rections for the future theoretical and experimental work 
in the field. On the whole, our results indicate that the 
effects of the Hubbard U are indeed crucial for the under- 
standing of magnetic properties of the CMR compounds. 
Qualitatively, this suggests that in addition to the fa- 
miliar double-exchange band-narrowing effects, the cor- 
related behaviour of spin- polarised large-U Hubbard car- 
riers should be recognised as an important mechanism 
underlying the physics of CMR compounds. 



II. SPIN WAVE THEORY 

Low-temperature spin-wave properties of double ex- 
change ferromagnets with U = are well understood (sec 
Ref. [27] and references therein). In parallel with Rcf. 
[27], we begin our treatment of the Hamiltonian, Eq. (1) 
with U ^ 0, with the standard Holstein-Primakoff trans- 
formation, followed by a canonical transformation of the 
form 
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where oL is the magnon creation operator, and N is the 
total number of lattice sites. The resulting Hamiltonian, 
Ti', takes form of a series in powers of 1/^/S <C 1, with 
the leading-order term, 
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(here, ^ means that the quasimomentum conservation 
law is obeyed). Since we will be interested in the leading- 
order (classical) spin wave properties, we will need only 
the two further terms in this series. Here, in addi- 
tion to the "usual" terms occurring already in the non- 
interacting model [27], 
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(where M^i.s stands for Wj^ ^ 
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In writing Eqs. (6-8), we omitted the terms contain- 
ing more than one spin-down fermion operator cj:. or cl 

(hence the "eff" above the equality signs). Due to the 
absence of spin-down electrons in the half-metallic ferro- 
magnetic ground state, these terms will not contribute to 
the quantities which are of interest to us here. 

We find it advantageous to choose the coefficients 
Wj: - in such a way that the leading-order single-particle 
electron-magnon scattering, Ti'j^, is cancelled by the aver- 
age contribution of Ti^j^ [see Fig. 1 (a)]: 
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Within the mean-field picture, condition (9) implies 
that the average number of magnons with a given mo- 
mentum p, 
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remains constant [28]. In other words, it represents the 
optimal choice in separating the two distinct branches of 
excitations (magnons and electrons/holes). 
Equation (9) is solved by 
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which at [/ ^ reduces to the familiar form [29], used 
earlier for the non-interacting case [27,29]. 

These expressions can be further simplified in the ex- 
perimentally relevant case of Js ^ ep, where ep = ji + td 
is the Fermi energy, measured from the bottom of the 
spin-up subband. In particular, we then find 
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If we also restrict ourselves to the case of small magnon 
momenta, p,p' <^ 1, the second term on the r. h. s. of 
Eq. (13) can be omitted, and we obtain 
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where v^ = dej: /dk is the electron velocity. 

The spin-wave energy, up^ is equal to magnon self- 
energy [30], which in turn can be evaluated perturba- 
tively (in l/\/S). In addition to the first-order contribu- 
tions from Ti'2 and ?i^2' there is a number of second-order 
corrections from Ti'i and Ti-'^i- Owing to the condition 
(9), these second-order terms cancel each other, with the 
sole exception shown diagrammatically in Fig. 1 (b). In 
drawing and evaluating this diagram, we make a drastic 
simplification of a mean-field type, corresponding to the 
Hartree decoupling of the interaction term in Eq. (4). 
Namely, we do not include any spin-up - spin-down elec- 
tron vertices and use the expression, 



Gi{uj,k) = l/{uj- ej:- Js + iJ. + iO- signer) (17) 

for the spin-down electron Green's function. The Green's 
function for a spin-up electron in the half-metallic case is 
given by the usual formula, G|(a;, k) — l/(a; — e^ + /i + 
iO ■ signcj) . 

The resultant expression for long-wavelength magnon 
dispersion takes the usual form, up — Dp^, where the 
spin-stiffness D is given by 
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The doping dependence of the spin-stiffness D for a two- 
dimensional system is illustrated in Fig. 2 (a). Here, 
the solid line shows our result, Eq. (18), for J = and 
experimentally relevant values Jn/t ~ 5, U/t = 16. The 
effect of Hubbard repulsion becomes clear from compari- 
son with the dashed line, corresponding to the Jn/t — 5, 
U = case [31]. We sec that in the presence of the on- 
site repulsion, the magnitude of D{x) increases, and the 
maximum is shifted towards x > 0.5. The dashed-dotted 
line corresponds to the Jn/t = 14.6, U = case, which 
is characterised by the same value of mean-field band- 
splitting Js as Jn/t = 5, U/t =16 system at the exper- 
imentally important value of electron density, x — 0.6. 
While numerically the two results at 0.5 < a; < 0.8 are 
relatively close, the dashed-dotted line possesses a larger 
slope, and still reaches maximum below quarter-filling, 
X = 0.5. For completeness, we note that the classical 
Jh -^ 00 result (dotted line) is symmetric and corre- 
sponds to the largest magnitude of D{x). We conclude 
that at a finite Jh, in addition to an overall increase in 
£), the inclusion of Hubbard repulsion leads to a relative 
increase of spin stiffness at x > 1/2, which is consistent 
with the experimental observation that the ferromagnetic 
tendencies in the CMR manganates are more pronounced 
in this doping region. We note that an earlier mean field 
study [7] of the effects of the Hubbard repulsion on the 
Curie temperature, Tc{x), suggested somewhat similar 
trends. 

The overall increase of spin stiffness originates from the 
mean- field effects discussed in Sect. I and corresponds to 
substitution Jh -^ Js in the last term of Eq.(18). Since 
Js > Jh, the (negative) pre- factor in front of the (pos- 
itive) integral decreases in comparison with the L'^ = 
case, resulting in the increase of D. At sufficiently low 
values of x, however, this tendency is counter-balanced 
by effects of the Hubbard correlations, which contribute 
the quantity U'^x{l — x) to the numerator of this pre- 
factor. The underlying physics will be discussed in Sect. 
HI below; here we merely note that as a result, the spin 
stiffness well below the quarter-filling, x < 0.2, may ac- 
tually be somewhat suppressed in comparison with the 
U = value. For the case of Jh — 5t, U = 16t, this takes 
place for x < 0.19 [see Fig. 2 (a)]. 



As for a quantitative comparison of our results for D{x) 
at x > 0.5 with the experimental data for the CMR man- 
ganates, this appears problematic due to a number of 
reasons: (i) available experimental results on the doping 
dependence of spin-stiffness [32] are still incomplete; (ii) 
it is known [27] that quantum corrections, not included in 
Eq. (18), lead to an appreciable rcnormalisation of spin 
stiffness magnitude; and (iii) the values of band (and or- 
bital) structure parameters and direct exchange integrals 
for particular compounds are known with a large degree 
of uncertainty. Nevertheless, it is adequate to say that 
qualitatively, both the overall profile and the magnitude 
of spin-stiffness, as given by Eq. (18), are consistent with 
the experimental results for D{x) within the metallic fer- 
romagnetic region, 0.2 < 1 — a; < 0.5. 

When spin stiffness [as shown in Fig. 2 (a) for J — 0] 
turns negative, the ferromagnetic ground state can only 
be stabilised by including a sufficiently strong direct fer- 
romagnetic exchange coupling, J < 0. On the contrary, 
a positive value of spin stiffness does not guarantee the 
stability of a uniform ferromagnetic ground state, since 
the latter might still be unstable with respect to phase 
separation (see Sect. IV below). 

Within the mean-field approach taken here, expression 
(18) is expected to hold at Js S> ep for all values of 
the ratio U/Jh, except for the Jh -^ case of a pure 
Hubbard model. In this case, the ionic spins are fully 
decoupled from the itinerant ones, and the leading-order 
(in 1/S) term in the magnon energy vanishes. Formally, 
at Jh -^ the second term on the r. h. s. of Eq. (13) 
cannot be omitted even for small p, and Eqs. (14-18) 
become invalid. We note that if the value of U is suf- 
ficiently large, the conduction electrons may still be in 
the fully spin-polarised ferromagnetic state as expected 
for a partially-filled Hubbard model below half-filling, 
a; < 1. While this is always the case within the present 
mean field treatment, actual identification of the stabil- 
ity region for a ferromagnetic state of a large-t/ Hubbard 
model remains an open problem [33,34]. Although this 
subject is well beyond the scope of the present work, it 
is important to note that (i) it is likely that over a broad 
range of doping values, the instability of the fully spin- 
polarised state of the Hubbard model {Jh — 0) results 
only in a partial reduction of magnetisation [33] , and (ii) 
it is possible that allowing for a small but finite value of 
Jh greatly enhances stability of the fully spin-polarised 
state (cf. Ref. [34]). 

In view of relatively small values of Jh, reported in the 
bandstructure calculations, it is important to study the 
crossover to the free-spin {Jh -^ 0) case in some detail. 
The leading-order (in both ep/ Js and 1/S) term in the 
magnon energy originates from 7^2 and in the absence of 
a direct coupling J has the form 
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Here, monicntuni p is allowed to take any value within 
the Brillouin zone. When Jh is sufficiently small, the 
second term in the denominator dominates [35] , provided 
that p^ ^ {JsJh)/{U\E\). The magnon energy then sat- 
urates at a constant value, cop = xJu / {2S), which is con- 
sistent with a physical picture of independent ionic spins 
Si subject to an effective magnetic field of the magnitude 
xJh/'^- The latter is created by the rigid ferromagnetic 
Fermi sea of the large-C/ Hubbard carriers [36]. This sit- 
uation is illustrated by a dashed-dottcd line in Fig. 2 
(h), corresponding to Jh /t = 0.2, U/t ~ 16 and x ~ 0.6. 
While for the experimentally relevant value of Jn/t = 5 
(solid line) spin-wave energy does not reach saturation, 
the effects of suppression of the magnon energy (in com- 
parison with the pure cosine Heisenberg law - see the 
upper dotted line, corresponding to Jh — > oo) are still 
felt near the zone boundary [37]. These become more 
pronounced (possibly leading even to a local minimum of 
spin-wave dispersion at the point p = {7r,7r}), if a direct 
antiferromagnetic coupling J > between the ionic spins 
is taken into account. The latter gives rise to an extra 
term, —2J{ej}+ td)/{tS), which should be added to the 
r. h. s. of Eq. (19). Thus, at the moderate values of 
Jh within the experimentally relevant range, the spin- 
wave dispersion in a double exchange ferromagnet with a 
large on-site repulsion U shows zone-boundary "soften- 
ing" in comparison with the nearest-neighbour Heisen- 
berg dispersion law. There has been an extensive theo- 
retical effort [38,39] directed at understanding this prop- 
erty, which is observed experimentally [40] in many (but 
not all [41]) CMR compounds. It is important that this 
generic feature is recovered within the present model, as 
suggested by earlier variational studies [39] of the effects 
of U on the magnon dispersion. 

We note that this zone-boundary softening effect oc- 
curs in both 2D and 3D, and is entirely due to the Hub- 
bard repulsion, U. Indeed, for any dimensionality d it can 
be shown [42] that a.t U — and x > 0.5, the magnon 
dispersion at sufficiently large Jh ~ e_F hardens towards 
the zone boundary. This is illustrated by the dashed line 
in Fig. 2 (b), representing the Jh = 5t, U ^ case [43]. 
One can see that at large momenta, the corresponding 
Heisenberg dispersion law, cop = 2D{td -\- ep)/t with the 
appropriate value of spin-stiffness D, indeed yields lower 
magnon energies (lower dotted line). 



III. EFFECTIVE ELECTRON-ELECTRON 
INTERACTION 

The original on-site Coulomb repulsion, as represented 
by the last term in Eq. (1), acts between electrons with 



anti-aligned spins. As the average number of spin-down 
electrons in a half-metallic double exchange ferromagnet 
at low temperatures is negligible, it might seem that the 
spin-up electrons remain non-interacting even in the pres- 
ence of the Hubbard U . However, as already discussed in 
Sect. I, at finite temperatures the on-site repulsion also 
gives rise to an interaction between electrons with the 
same sign of spin projection. The presence of this novel 
interaction, Vg//, is clear, e.g., from the form of the oper- 
ator 7i^2 [see Eq. (8)], which at finite T can be averaged 
over the equilibrium magnon distribution. Another con- 
tribution [44] to Veff originates from the second-order 
processes, involving various combinations of terms from 
Ti'i and Ti'^n Eqs. (5) and (7). As in Sect. II above, the 
condition (9) leads to a massive cancellation among these 
second-order diagrams [cf. Fig. 1 (a)] , with the only two 
surviving terms shown in Fig. 1 (c). 

In Fig. 1 (c), vertices correspond to the magnon- 
electron interaction, Ti-'n, and the solid lines are finite- 
temperature electron Green's functions, 

(20) 

where we again used the Hartree mean-field form for Q^ 
[cf. Eq. (17)]. At low temperatures, T <$^ D, only the 
long- wavelength magnons are present, and the magnon 
Green's function can be written as Gm^ ~ iC, — Dp^. Here, 
the spin stiffness D can be evaluated with the help of Eq. 
(18) or taken directly from the low-temperature neutron 
scattering measurements. Furthermore, at low tempera- 
tures Eqs. (15-16) may be used in place of Eqs. (7-8). 

To leading order in ep / Js a-nd l/S*, the net result for 
the vertex function, Fjj 43, of two spin- up electrons scat- 
tering is given by the expression 
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which should be anti-symmetrised [45] with respect to 
the velocities and frequencies of the outgoing (incoming) 
electrons, wi^2 = dep^^/dpi^2 and ^1,2 (^^3,4 and C3,4)- In 
Eq. (21), the first term in brackets represents the first- 
order contribution of the operator 7i^2' ^'l- (16), whereas 
the other two terms come from the two diagrams shown 
in Fig. 1 (c); Np — [exp(Dp^/r) — 1]^^ is the average 
magnon occupation number, Eq. (11). 

The retardation (frequency-dependence) effects in V^^ 
become noticeable only on a very large electron energy 
(frequency) scale [46] of \C,i\ ~ Js, and can be omit- 
ted whenever only electrons with energies near the Fermi 



level are considered. In this case, the effect of electron- 
electron scattering as described by the vertex F^^^ is 
equivalent to that of an effective electron- electron inter- 
action of the form 



with 



Veff = 



U[Js + U{2x - 1)] 
SdSJlN 
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where we also used the fact that the long-wavelength 
magnon dispersion is isotropic. In the case of small 
U ^ Jht^f, our mean-field results, Eqs. (21) and (22), 
can be re-derived within the perturbation theory in U. 
As expected on physical grounds, Veff vanishes also in 
the U ^ 00 limit, when the double occupancy on-site is 
forbidden. The momentum integral occurring on the r. 
h. s. of Eq. (22) can be easily evaluated. 
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in 2D, 



where C(5/2) ~ 1.34 is the Riemann's zeto- function. We 
note that the quantity (23) can also be expressed macro- 
scopically as the thermal average of 
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where M" are the three components of local magnetisa- 
tion [47], M. This shows that the appearance of Veff 
is indeed a direct consequence of the misalignment of 
neighbouring spins (which in turn is due to the thermal 
fluctuations; cf. Sect. I). 

Although the precise form of Veff in Eq. (22) obvi- 
ously has only a mean-field validity, we emphasise that 
qualitatively this effect, which has a clear physical ori- 
gin, will survive in an exact treatment. If anything, the 
mean- field approach yields a smaller magnitude of Veff : 
indeed, Eqs. (20) overestimate the energy of a spin-down 
electron, which enters into the denominators of diagram- 
matic expressions in Fig. 1 (c). 

Thus, we conclude that with increasing temperature 
T, there arises an effective interaction between the spin- 
polarised carriers in a double exchange ferroniagnet with 
U > 0. For the purposes of order-of-magnitude es- 
timates, one can assume that electron dispersion is 
isotropic, H w p/m^, where m^ is an effective mass of 
electron (or hole). The interaction, Eq. (22), then takes 
form of a simple p-wave scattering. 
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In real space, the effective interaction, Eq. (22), takes 
form 
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Here, for each lattice site i a summation over its 2d near- 
est neighbours (labeled i + A) is performed. The first 
term in Eq.(27) contains a product of carrier densities; 
it renormalises the "bare" repulsion between electrons on 
the neighbouring sites, which is present in reality but not 
included in our model, Eq. (1). In the case of the CMR 
manganates, this bare repulsion is expected [48] to be of 

the order of Vnn ~ 0.05 eV; on the other hand, the value 
of tV(T) as found from Eq. (26) with U = 8cV,t = 0.5 
eV, Jh = 2.5 eV and x = 0.6 is t^V{T) ~ IQ-'^iT/Tc)^/^ 
(in units of eV, for the 3D case; we also assumed that 
the Curie temperature Tc is of the order of the zero- 
temperature spin-stiffness, D). We thus see that in the 
manganates, the contribution of Veff, Eq. (27), to the 
nearest-neighbour repulsion remains negligible even for 

The situation is likely to be very different for non- 
manganate magnetic semiconductors exhibiting CMR. It 
is expected [49,50] that in the lightly doped {x ^ 10^^) 
EuSe , EuTe, EuO or undoped (semimetallic) EuBg the 
parameters of the Hamiltonian, Eq. (1) can be roughly 
estimated as C/ = 7 eV, Jh = OAeV and t = 0.5 eV. Tak- 
ing also into account the magnitude of Eu spin (7/2), we 

find t^V{T) 0.2(r/Tc)^/2 in units of cV. Thus, the 

effects of temperature-dependent renormalisation of Vnn 
in these compounds may be appreciable. 

We note that in the latter example, the sign of V{T) 
is negative, corresponding to an effective attraction be- 
tween electrons on the neighbouring sites. This clearly 
contradicts the simple physical picture outlined in Sect. 
I (based on the small-[/ perturbative considerations). In- 
deed, as can be seen from Eq. (26), at sufficiently large 
values of U the effective interaction can in fact become 
attractive. This occurs when the second term in the nu- 
merator, U{2x — 1), which is negative for x < 0.5, dom- 
inates over the first one. We suggest that this second 
term, which originates from the two diagrams shown in 
Fig. 1 (c), is due to the many-body effects. 

Indeed, the large-C/ partially-filled Hubbard conduc- 
tion band has (at least within the mean-field picture - see 
the discussion in Sect. II above) the ferromagnetic ten- 
dencies of its own, which are unrelated to the ionic spins 
and double exchange. In the Jh = case, the ferromag- 
netic, fully spin-polarised state of the carriers (uncoupled 



to the ionic spins) corresponds to the largest bandwidth 
and hence to the lowest kinetic energy. Let us now as- 
sume that the value of Jh is finite and the ionic spins are 
fixed in a certain configuration which is not perfectly fer- 
romagnetic (in the present context, the deviation of the 
ionic spin configuration from the ferromagnetic ground 
state is due to the thermal fluctuations). The electron 
bandstructure is then determined by a competition be- 
tween the Hubbard band ferromagnetism (which favours 
a uniform ferromagnetic alignment of the carrier spins 
and hence decoupling from the ionic spin background) 
and the Hund's rule coupling (which tends to align the 
carrier spins locally with the ionic ones, leading to the 
double exchange band narrowing) [51]. Not surprisingly 
for an interacting many-body system, the resulting band- 
structure (along with the local carrier spin direction and 
the overall strength of band ferromagnetism) therefore 
depends on the bandfilling x. As a result of electron- 
electron interaction, the mean-field bandwidth of spin- 
polarised carriers (somewhat suppressed due to the dou- 
ble exchange mechanism) then increases whenever the 
bandfilling of electrons (or holes) is increased, reaching 
the maximum at x = 0.5 (where the effects of Hubbard 
band ferromagnetism arc most pronounced) . The kinetic 
energy of an electron at a sufficiently small value of x may 
thus be lowered if another electron is present nearby (a 
local increase of the carrier density), leading to the effec- 
tive attraction. 

This effect is likely to have important physical conse- 
quences in the case of Eu-based magnetic semiconduc- 
tors, as the effective attraction will further improve the 
stability of microscopic droplet-like areas with increased 
carrier density and enhanced ferromagnetic order (fer- 
rons [2]), which were argued [2,52] to play a key role in 
these compounds. If the bare nearest- neighbour Coulomb 
repulsion is not too strong (which is expected), taking 
the effective interaction into account may in fact lead 
to the total nearest neighbour interaction being attrac- 
tive above a certain temperature. This in turn might 
signal an instability of the homogeneous ferromagnetic 
state and possibly the formation of ferrous. This ques- 
tion obviously calls for further investigation. 

The second term in Eq. (27) is a variant of corre- 
lated electron hopping. While the effects of this term in 
the present context are not immediately clear, we note 
that correlated hopping represents a much-studied ex- 
tension of the Hubbard model. In fact, an effective in- 
teraction somewhat similar to (27) (although involving 
fermions with antiparallel spins) has been obtained in the 
past [53] within a mean-field approach to a \a,vge-U Hub- 
bard model. In the case of CMR manganates, where the 
quantity ^(r) is very small, the effects of the correlated 
hopping term in Eq. (27) are expected to be negligible. 
Again, the situation can be very different for the CMR 
magnetic semiconductors with small carrier densities. 

We note that while the result (27) apphes only for de- 



generate magnetic semiconductors with finite carrier den- 
sities at T ^ 0, we expect a similar magnon-mediated 
effective electron-electron interaction to arise in the non- 
degenerate (undoped) case as well. 



IV. PHASE SEPARATION AND PHASE 
DIAGRAM AT T=0 

Although spin stiffness D{x) and spin- wave energy, Wp-, 
discussed in Sect. II, are important and much-studied 
quantities characterising magnetic properties of double 
exchange ferromagnets, their behaviour is not expected 
to yield any conclusive information on the stability of 
the ferromagnetic state at low temperatures. Indeed, for 
large Ju and C/ = it is possible to verify (see, e.g., Ref. 
[18] and references therein) that, in a marked difference 
from conventional isotropic insulating magnets with com- 
peting interactions, the zero-temperature magnon spec- 
trum of a double exchange magnet in the homogeneous 
ferromagnetic state does not soften when the latter is 
rendered unstable due to a change in the balance be- 
tween ferro- and antiferromagnetic tendencies in the sys- 
tem. This is because phase separation, which is a generic 
phenomenon found in both experimental and theoreti- 
cal studies of the CMR manganates and doped magnetic 
semiconductors [2,17] always preempts a second-order, 
spin wave-mediated phase transition. On the theory side, 
there is little doubt that this situation persists in the 
finite- J//, U > Q case in two and three dimensions (in 
particular, this is known to be true in the Jh —^ 0, large- 
U case of the pure Hubbard model [54]). In the present 
section, we will see how the presence of on-site Coulomb 
repulsion U affects the physics of phase separation and, 
in particular, the low-temperature phase diagram. We 
note that the effects of Coulomb repulsion on phase sep- 
aration in this model were studied in Ref. [4] within a 
somewhat different mean-field approach. 

At a given value of electron band filling x [and the cor- 
responding chemical potential, iJ.{x)], the homogeneous 
ferromagnetic phase of a double exchange magnet is un- 
stable with respect to phase separation whenever its ther- 
modynamic potential. 
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is larger than the thermodynamic potential Q,p of an- 
other homogeneous phase P, calculated at the same value 
of /i. The total energy of the system can then be lowered 
if areas of this new phase are formed within the bulk 
of the ferromagnct [55]. Therefore in order to find the 
stability region of the homogeneous ferromagnetic phase, 
one has to identify the relevant phases and evaluate the 
dependence of their respective thermodynamic potentials 
on the parameters of the Hamiltonian, Eq. (1). 



For any given values of Jh and U , it is convenient to 
characterise the stabihty of the homogeneous ferromag- 
netic state by the critical value of direct antiferromag- 
netic exchange coupling, Jcr{x), above which the system 
becomes phase-separated. The calculation thus proceeds 
as follows: the equation, 



ftpAiiKx), Jp] = flp[n{x), Jp] , 



(29) 



is solved for several possible phases P, yielding the cor- 
responding values of J = Jp{x,U,Jh)- At a fixed x, 
the value of Jcr is then given by the lowest Jp. Such a 
procedure clearly has a variational validity as it does not 
imply an existence of a rigorous proof that Jcr cannot be 
lowered further by broadening our selection of phases P. 
On the other hand, the condition J > Jcr is obviously 
sufficient for the phase separation to occur. We shall 
now turn to the two possible antiferromagnetic phases 
considered in the present paper. 

Based on the numerous results for the U = case 
[2,18,17], one expects that a phase separation into ferro- 
magnetic and the usual Neel antiferromagnetic (G-type 
antiferromagnetic) phase, corresponding to the wave vec- 
tor of {tt, tt, tt} (in 2D, {tt, tt}) takes place in the vicinity 
of endpoints, x = and x = 1. To leading order in 1/S, 
one can use the classical formalism {see, e.g., Eqs. (9- 
10) of Ref. [15]} to re-write the Hamiltonian in terms of 
the new fermions di^- (and di^), whose spins are aligned 
(antialigned) with the local ionic spins Si in the G-type 
antiferromagnet : 

t 
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Here, a standard Hartree mean-field decoupling has been 
carried out in the last term, with Xa = {dladia) de- 
noting the average spin-up and -down fermion densi- 
ties. Upon Fourier transformation and diagonalisation 
(see Appendix A), Eq. (30) takes form 
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Here, /I are the new fermion creation operators, the 

momentum summation is performed over the full (fer- 
romagnetic) Brillouin zone, and the carrier spectrum is 
given by 



^h'ik) = ST\l^(jPr + {ej:r, 



(32) 
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These equations are valid also at [/ = 0, in which case 

(G) 

Jg = Jh and 5 = 0. For general values of U, the 

(G) 

quantity Jg represents the Stoner-type mean field band 
splitting in the G-antiferromagnctic phase, whereas 6 
merely renormalises the chemical potential, ^—-^Jh- The 
closed system of mean field equations for these parame- 
ters of the G-antifcrromagnetic phase at a given band 
filling, xq = x^ -\- x^, is presented and discussed in Ap- 
pendix A. As we are interested only in finding the phase- 
separation instability of the homogeneous ferromagnetic 
state, we do not need to solve these mean field equations 
for general values of xg ■ The latter will rather be deter- 
mined by ji, which in turn is related in a usual way to 
the band filling x of the homogeneous double exchange 
ferromagnet. 

When the value of x is small, the chemical potential 
in the ferromagnetic phase lies below the bottom of the 
lower band of G-antiferromagnet [given by Eq. (32) with 
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and the thermodynamic potential of the Neel phase, Qg, 
equals —dJ. As the value of x increases, the inequality 
(34) is eventually violated, giving rise to a non-zero car- 
rier density xg (with X| > x^) in the antiferromagnetic 
phase. 

For the realistic values of parameters in the 2D case, we 
find that the G-antiferromagnetic phase with partially- 
filled band, > xg > 1, is not relevant in the context 
of phase separation. This is because within the corre- 
sponding range of values of x [and /i(a;)], the critical 
value of superexchange coupling, Jg{x), is larger than 
the one which corresponds to the second-order spin- wave 
transition, Jaw{x) ~ Dq{x)S (where Dq is the spin stiff- 
ness evaluated at J = 0). The peculiar properties of the 
partially-filled case, as discussed in Appendix A, are in- 
teresting on their own, and may also prove relevant in 
another context. Here, we will turn to the case when the 
chemical potential ix{x) lies above the top of the filled 
G-antifcrromagnctic band. 
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resulting in xq = 1. Given the reahstic (large) magnitude 
of [/, the value of Jg. « Jh -\-U \s then also large, and 
the mean field equations are readily solved with the help 
of the t/ Jg -expansion [which was used also in writing 
Eq. (35)]. We find an expression 
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which can be used to determine the critical value of su- 
perexchange Jg(x) whenever the phase separation into 



the partially filled ferromagnetic and filled {xq = 1) Neel 
states becomes possible. As we will see below, this hap- 
pens only well above the quarter-filling of the ferromag- 
netic band, x > 0.5, where the inequality (35) is clearly 
satisfied. 

Another phase which seems to be ubiquitous in the 2D 
case is the A-type antiferromagnetic one, characterised 
by the wave vector {7r,0}. The mean field theory of 
the A-antiferromagnetic phase is formulated along the 
same lines as for the Neel antiferromagnet above, where 
Eq.(32) is now replaced with 
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)-tcosk^T JjiJs^^)^ + ^^ cos2 ky , (37) 

and the (negative) superexchange term in the Hamilto- 
nian is cancelled. In spite of the higher superexchange 
energy, the ^-type phase becomes more profitable than 
the Neel one because it allows for a larger gain in the 
kinetic energy of carriers. Thus, it is clear that the A- 
phase becomes relevant only when the carrier density in 
the A-type ferromagnet ^^[/^(a;)] (and similarly, the hole 
density, 1 — x^) is not too small. In case of large U, this 

(A) 

means that the value of Stoner bandsplitting, Jg , is 
much larger than the hopping matrix element, t. 

Therefore while for moderate values of U the full sys- 
tem of mean field equations for Jg and xa must be 
solved, and then the thermodynamic potential. 
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[where n- is the Fermi distribution function, corre- 
sponding to the dispersion law (37) and chemical po- 
tential /i(x) — Jh/'2], must be evaluated numerically, in 
the large-C/ case it is sufficient to retain the leading-order 
terms in t/ji . In this way we obtain for U ':$> t, 
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Eqs. (39-40) hold for \iJ.(x)\ < t, which at large U corre- 
sponds to partial filling of the A-antiferromagnetic band. 
The zero-temperature phase diagram for a 2D system 
with Jn/t ~ 5 and U ~ is shown in Fig. 3 (a). 
Here, the (upper) dotted line corresponds to the long- 
wavelength spin- wave softening, Jsw{x) = Do{x)S; at 
X < 0.31, it is preempted by spin wave instability at 
p = {tt, tt} (lower dotted line) . Phase separation instabil- 
ities involving Neel and A-type antiferromagnetic phases 
are represented by solid and dashed lines, respectively. 



We see that the stability region of the uniform ferro- 
magnetic phase is heavily shifted towards the electron- 
doped end, X < 0.5; at the experimentally relevant val- 
ues of J/t ~ 0.015 {J/t ^ 0.02), the homogeneous fer- 
romagnetic state is unstable everywhere at 1 — a; < 0.34 
(1 — a; < 0.55). As explained in the Introduction, it is at 
X > 0.5 that the broad low-temperature ferromagnetic 
region is found for the CMR manganates, and such an 
instability clearly contradicts this experimental observa- 
tion. 

The situation changes in the presence oi U — 16t, a,s 
shown in Fig. 3 (b). Here, the dotted line represents 
our result for spin stiffness, Eq. (18); due to the zone- 
boundary softening effect discussed in Sect. II, it is ex- 
pected that for all carrier concentrations, the spin- wave 
instability at p = {tt, tt} corresponds to a slightly lower 
value of superexchange. The latter is not shown in Fig. 
3 (b) because our result for the magnon spectrum, Eq. 
(19), does not include the subleading (in t/Js) term and 
therefore should not be compared with other quantities 
plotted here. The Neel phase separation instability (solid 
line) is plotted using Eqs. (34), (A12), and (36). In or- 
der to improve accuracy for /x(x) < — 1 (corresponding to 
X < 0.19), we solved the full system of mean field equa- 
tions for the A-type phase. Nevertheless, we note that 
Eq. (40) yields accurate results in the region where the 
phase separation into ferromagnetic and A-type antifer- 
romagnetic phases is possible. 

We see that for U/t = 16 and J/t = 0.015 {J/t = 
0.02), the region where ferromagnetic phase is unstable 
is shifted to 1 — a; < 0.21 (1 — a; < 0.23), so that a sub- 
stantial stability region is now left for the ferromagnetic 
phase at X > 0.5, in agreement with experimental results. 
Admittedly, the presence of broad stability region of fer- 
romagnetic phase at a: < 0.5 is at variance with experi- 
ments and indicates a deficiency of either our variational 
procedure (i.e., our choice of possible phases is too nar- 
row), or our simplified model, Eq. (1) [56]. Indeed, in 
the 3D case the experimental data [10,57] show that the 
phase diagrams for different perovskitc compounds dif- 
fer in the electron-doped {x < 0.5) region, suggesting 
tne sensitivity to details of crystalline surrounding and 
perhaps the importance of orbital structure. 
(40) The inclusion of Hubbard interaction in principle could 
have brought about new charge-ordered antiferromag- 
netic (or ferrimagnetic) phases, which would not occur 
at C/ = 0. While we tried to look into this possibil- 
ity, we could not identify any such phases that would be 
stable within the experimentally relevant region of pa- 
rameter values. This difficulty was encountered also by 
other workers in the field [4], who could stabilise charge 
ordering only upon including a large intcrsitc Coulomb 
repulsion. 

It is however worth mentioning that some of the more 
complicated phases which are relevant for phase separa- 
tion at C/ = do show charge ordering owing to inequiv- 
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alence of different lattice sites. This is exemplified by the 
chain phase [18], shown schematically in Fig. 4. In Fig. 
3, the critical value of superexchange, corresponding to 
phase separation into ferromagnetic and chain phases, is 
shown by the dashed-dotted line; for the U = case it 
was calculated exactly, whereas for U = 16t wc used the 
large Stoner bandsplitting expansion similar to the one 
described above for the y4-antifcrromagnetic phase. We 
see that when U is included, the phase separation into 
ferromagnetic and chain phases in the hole-doped region, 
X > 0.5, becomes impossible. This is in line with the gen- 
eral expectation that Hubbard repulsion disfavours charge 
ordering at x > 0.5, when there is less space between 
electrons. 

The present mean-field treatment allowed us to arrive 
at important conclusions regarding the effects of Hub- 
bard repulsion on the possible instabilities of the homoge- 
neous ferromagnetic phase. Nevertheless, we stress that 
the full zero-temperature phase diagrams for the model 
(1) both in 2D and in a much more cumbersome 3D case 
are still lacking, and should only come from numerical 
experiments. This is a challenging problem, as some of 
the phases involved can be expected to have relatively 
large unit cells. 



V. CARRIER DENSITY OF STATES NEAR THE 
FERMI LEVEL 

While the physical nature of carrier transport and 
magnetotransport in the CMR compounds near the Curie 
temperature remains largely unknown, it may be possi- 
ble to single out an equilibrium property which is most 
closely related to the CMR phenomenon. It appears that 
such a property is a broad depletion of carrier density 
of states near the Fermi level as observed in photoe- 
mission/absorption in the CMR manganates [25]. This 
decrease of the density of states, visible already deep 
in the ferromagnetic phase, becomes progressively more 
pronounced as the temperature approaches Tq, around 
which there is no spectral weight left at the Fermi level 
within the accuracy of the experiments. These results 
were subsequently confirmed by the tunnelling measure- 
ments [26] ; it was suggested that the hard gap opening at 
T w Tc is responsible for the peak of the resistivity. The 
relevance of these gap or pseudogap phenomena for trans- 
port is underlined by the fact [26] that the "transport 
gap" seen in the activated temperature dependence of 
resistivity at T > Tc is roughly of the same order of mag- 
nitude (tenths of eV) as the width of the density of states 
depletion. Furthermore, it seems possible that in the case 
of the non-manganate lightly doped or non-degenerate 
magnetic semiconductors, the well-known giant red shift 
of the optical absorption edge [58] (with optical gap de- 
creasing as the temperature is lowered through Tc] see 
the discussion in Ref . [2] ) may act as a counterpart of the 



temperature-dependent pseudogap observed in the man- 
ganates. 

For the case of the CMR manganates, high-resolution 
tunnelling measurements have recently been extended 
[24] down to liquid helium temperatures, revealing a no- 
ticeable, albeit narrow, depletion of the density of states 
near the Fermi level at T = 4.2K. Although the con- 
clusive evidence is still lacking, it is most reasonable to 
expect that it is this feature which with increasing tem- 
perature evolves into the broad pseudogap observed near 
T = Tc- It is therefore important that a proper descrip- 
tion of the low-temperature ferromagnetic state of the 
CMR compounds should include this anomaly. 

Within the mean-field picture advanced in the present 
paper, at low temperatures and within the relevant dop- 
ing range of 1 — a; « 0.3, the system is assumed to be 
in a homogeneous, half-metallic ferromagnetic state; de- 
viations from ferromagnetic ordering (spin waves) freeze 
out at T <C Tc- The resistivity of the sample is then 
determined by impurity scattering, and it is only a com- 
bination of the latter with the usual Coulomb repulsion 
(including the long-range component), which could result 
in any density of states feature near the Fermi level. In- 
deed, it is known from the work of Altshuler and Aronov 
[59,60] that electron-electron interactions in a diffusive 
conductor generate an anomaly in the tunnelling density 
of states, centered on the Fermi energy. We note that 
since the magnitude of effective interaction Veff{T) in 
the manganates is much smaller than the corresponding 
(nearest-neighbour) term Vnn in the Coulomb repulsion 
(see Sect. Ill), the effects of Veff(T) on the density of 
states can be omitted altogether. 

We use the standard expression for the change in the 
density of states near the Fermi level [59,60], 



dv^ 



1 



M 



1/2 



V27r2 ihD,r/2 ' 



(41) 



where the overall pre-factor has been multiplied by two 
in order to account for half-metallicity of the system [61]. 
Assuming t ~ 0.5eV and using the values p ^ 161/ir2cm 
for the resistivity of Lao.rCao.sMnOs [24] and a ^ 4A 
for intersite distance, we estimate the diffusion constant 
De as De ^ [pe'^v{eF)r^ ^ 6aH/{pe'^) - 7cmVs. The 
resulting estimate. 



Si^ie) 



0.03 



Ml 



t 



(42) 



is an order of magnitude smaller than the experimental 
results of Ref. [24], which show a 15% change in Si^/i^ for 
|e — /i| ^ 0.075eV. Furthermore, based on Eq. (41) one 
expects that the relative change in the density of states 
for Lao.75Sro.25Mn03 (characterised by smaller values of 
resistivity and by a larger bandwidth) should be about 10 
times less than in the case of Lao.yCao.sMnOs, whereas 
experimentally the two curves differ by a factor of the 
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order of 2. In addition, the experiments [24] yield bvjv ex 
(e — /i)^ within a relatively broad energy range of |e — 
p\ < 0.02cV; this is in contrast with the standard theory 
[59,60], which predicts a crossover to the square root law 
[cf. Eq. (41)] at |e - ^| ~ T. Thus, we find that the 
Altshuler - Aronov mechanism cannot possibly account 
for the Fermi-level density of states depletion in the CMR 
manganates even at low temperatures. 

We therefore conclude that the low-temperature ferro- 
magnetic state of the CMR manganates is characterised 
by strong electron correlation effects. Since these are 
certainly not captured within the present mean-field ap- 
proach this does not necessarily signify the deficiency of 
our simplified model, Eq. (1). While we plan to in- 
vestigate this question in more detail in the future, we 
emphasise that our conclusion on the correlated nature 
of the low-temperature ferromagnetic phase of the CMR 
manganates is likely to be model-independent. In other 
words, an adequate generic model of the manganates, 
whether or not it involves orbital, lattice, etc. degrees of 
freedom, must necessarily take proper non-perturbative 
account of electron-electron interactions. 



VI. CONCLUSION 

In the present paper, we were concerned with the ef- 
fects of the strong on-site Coulomb repulsion which is 
present in the CMR compounds but often overlooked in 
the theoretical treatments. By treating the model within 
the mean field approach, we were able to resolve some 
apparent discrepancies between the generically observed 
low-temperature properties of these compounds and the- 
oretical results for the (non-interacting) double exchange 
model at the appropriate values of Hund's rule coupling 
strength. These properties include the doping depen- 
dence of spin stiffness, and the "zone-boundary soften- 
ing" of magnon spectrum (Sect. II) which has attracted 
much attention from both theorists and experimentalists. 
The underlying physical mechanisms are two- fold and in- 
clude both the interaction-induced increase in the effec- 
tive band splitting (Sect. I) and the correlated physics 
of the strongly interacting Hubbard carriers (Sect. III). 
In addition, we showed that a novel, magnon-mediated 
effective electron-electron interaction arises in these sys- 
tems at finite temperatures (Sect. III). While for the 
CMR manganates the strength of this interaction re- 
mains negligible, it is expected that it is much more im- 
portant in the case of Eu-based magnetic semiconductors 
exhibiting CMR. 

By considering the stability of the ferromagnetic state 
against phase separation, we were able to show (Sect. 
IV) that inclusion of the Hubbard repulsion alleviates an- 
other disagreement between the theory and experiment, 
resulting in a sizable stability region of the ferromagnetic 
state above half electron filling, x > 0.5. Regarding the 



phase diagram, the question of identification of the rele- 
vant phases and finding the precise domain of the ferro- 
magnetic phase (especially in the electron-doped region, 
X < 0.5) remains open and calls for further theoretical 
investigations, in particular numerical ones. At the same 
time we note that the underlying physics consists in a 
competition between many phases with very close values 
of thermodynamic potential, and the outcome is guar- 
anteed to be strongly dependent on the details of band 
structure, lattice/orbital properties and interactions in 
a particular compound. Therefore, while understanding 
the details and implications of phase separation in the 
CMR compounds (including both thermodynamic and 
transport properties) presents a broad and fascinating 
problem, it is not obvious that these details are directly 
related to the generic features of the CMR phenomenon 
itself. 

The satisfactory results of our mean field approach, as 
sketched above, all have to do with the integral quanti- 
ties, involving summation over the entire Fermi sea. It is 
precisely this effective averaging that makes our Hartree- 
Fock decoupling scheme a relatively reliable tool in this 
case. The situation changes drastically when this ap- 
proach is used to address other issues, such as the be- 
haviour of the carrier density of states near the Fermi 
level (Sect. V). In this case, the usage of the Hartree- 
Fock approximation (which yields the effective Stoner 
band splitting Js much larger than the Fermi energy, 
leading to an assumption that the system is half-metallic 
with no minority-spin carriers) is the probable cause of 
our inability to reproduce the experimentally measured 
[24] depletion of the density of states. Assuming that the 
results reported in Ref. [24] are sufficiently generic, this 
failure may have far-reaching conceptual consequences. 

As mentioned in Sect. V, it is expected that the low- 
temperature depletion of the density of states at the 
Fermi level is a temperature-dependent feature, which 
with increasing temperature evolves into the pseudogap; 
this temperature dependence is in turn expected to cru- 
cially affect transport properties of the system both near 
Tc and at low temperatures. The failure to recover 
the low-T Fermi-level feature in the density of states 
within a theoretical treatment based on the picture of 
a half-metallic homogeneous ferromagnetic phase (Sect. 
V) should lead to questioning the experimental relevance 
of the many available calculations of the low-temperature 
resistivity in double exchange ferromagnets, which are 
based on similar assumptions. 

On the other hand, one should not overlook the evi- 
dence, coming both from band structure calculations [14] 
and the experimental observations [62], which points to 
the presence of carriers in the minority spin subband in 
the CMR manganates even at low temperatures. The- 
oretically, this may be possible due to relatively small 
values of Hund's rule coupling; the spin-down (minor- 
ity) electrons [63], if present, will form polaron-like lo- 
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calised states, accompanied by a reduction in the spin-up 
electron density within the area of the polaron (this re- 
duction will in turn reduce the Coulomb/Hubbard inter- 
action energy). These localised spin-down carriers will 
clearly lead to an enhanced spin-up electron scattering 
and will also affect the density of states at the Fermi 
level [64]. As explained in the Introduction, these cor- 
related effects are beyond the mean-field approach used 
in the present paper; in our view, this scenario definitely 
merits further attention, especially as there is now some 
recognition [65] that other avenues of theoretical inves- 
tigation of the CMR and related phenomena may have 
proved unpromising. 

On the experimental side, we suggest that some key 
measurements still have to be performed in order to clar- 
ify the issues under discussion here. These fall into three 
categories: 

(i) Detailed investigations of the energy dependence of 
the density of states as a function of temperature and 
magnetic field. This includes tunnelling and optical mea- 
surements for various chemical composition and dop- 
ing levels, within the entire temperature range, and 
would help to clarify the relationship between the low- 
temperature density of state depletion of the Fermi level 
[24] and the pseudogap observed in the near-critical re- 
gion [25,26], as well as confirm the relevance of these 
phenomena for transport and magnctotransport. 
(ii) Systematic investigation of the interplay between 
pseudogap and other phenomena, in particular those re- 
lated to the unusual spin correlations found in the CMR 
manganates, such as the central peak observed in the 
inelastic neutron scattering [66] and critical behaviour 
of spin-stiffness [67]. In particular, it can be expected 
that understanding the nature of spin dynamics at el- 
evated temperatures would shed light on the structure 
of electron states involved in the pseudogap formation. 
The relationship between pseudogap and phase separa- 
tion should be clarified as well. 

(iii) Lastly, what appears an important theoretical and 
experimental problem is to identify the "common denom- 
inator" between the structure and properties of the CMR 
manganates [1] and those of the CMR magnetic semicon- 
ductors (Eu-based) and spinels [2]. Understanding what 
these CMR compounds have in common may be of much 
help in constructing a minimal theoretical model for the 
CMR compounds. In connexion to this, we recall that the 
familiar argument concerning the relative unimportance 
of U in systems with low carrier densities rests upon an 
assumption that carrier distribution throughout the sam- 
ple is uniform on the microscopic scale. The latter is not 
expected to be true, both for the CMR manganates and 
for ferromagnetic semiconductors at elevated tempera- 
tures, T ~ Tc- It is therefore possible that the effects 
of U, discussed in this paper, are important in both sys- 
tems. 
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APPENDIX A: MEAN FIELD EQUATIONS FOR 

THE G-TYPE ANTIFERROMAGNETIC PHASE 

AT T=0 

Here we outline the necessary details of the Hartree 
mean field scheme, as applied in Sect. IV to the G- 
antiferromagnetic phase. Throughout this appendix, we 
use units in which the hopping constant, t, is equal to 
unity. Eq. (30) is diagonalised by 
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with Vii — \{Js Y + e|. Assuming that the system is 
still haft-metallic, i.e. that the chemical potential lies be- 
low the bottom of the spin-down antiferromagnetic band 
[cf. Eq. (32)], we find 



=Ta 



1 1 ^ 

2 ^ 2 ' 



(G) 
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(A2) 



where ni. = (./I .fj:^) is the appropriate Fermi distri- 
bution function. Together with Eqs.(32- 33), Eq.(A2) 
forms a closed system of mean field equations for a ho- 
mogeneous G-antiferromagnetic phase at fixed xq- 

We note that the Fermi surface in a partially-filled 
spin-up band has two sheets, corresponding to different 
signs of er. Thus, the quantity 



{4^d^l) = -T7^Y. 



(G) 

^£"fc- 

^fc- 



is always equal to zero, showing the consistency of mean 
field decoupling used in Eq. (30). 

When the carrier density in the ferromagnetic phase 
becomes sufficiently large for the inequality (34) to be 
violated, there arises a non-zero carrier density, xg — 
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Xi + X| in the Neel phase. This in turn leads to an in- 

(G) 

crease of quantities S and Jg in Eq. (32), resulting 
in the upward shift and narrowing of the antiferromag- 
netic band. This effect is more pronounced when the 
value of U is sufficiently large, in which case the energy 
difference between chemical potential and the bottom of 
G-antiferromagnetic band, 



C-M(2:)-^Jif-<5+\/^(J^)2 + rf^ 



li^P' 



(A3) 



and the band-filling. 
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in 3D, 



in 2D, 



(A4) 



remain small within an extended range of values of x. In 
this large-C/ limit, it is possible to solve the mean-field 
equations analytically. 

When the band filling is small, xg <C 1, Eqs. (A2) and 
(33) yield 



Jf^ = Jh + Ujf^ 



or 



j(G) _ Jh ^ ^ 
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(A6) 



Substituting this into Eq. (A3) [and also using Eq. (33) 
for S] , we obtain 



C+ljH-Kx) = J\j!i+d'(^l-luC^ (A7) 

We are interested in the large-t/ situation when C. is small 
and can be omitted on the 1. h. s. [see below, Eq. (A13)], 
in which case we find, to leading order. 
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With the help of Eqs. (A6) this in turn yields 
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While at /x(x) 

a further increase of x 



Ho, Eq. (A9) yields J^*^^ = Jh, with 



[and hence /i(x)] towards the 

(G) 



quarter-filling, iJ,{x) = 0, the value of Jg increases, and 

the bandwidth [of the order of cP/{jf^f , see Eq. (32)] 
decreases, with the effect that the bottom of spin- up band 
in the G-antiferromagnetic phase remains pinned imme- 
diately below the chemical potential /i(a:;) — \Jh- As a 
result, band-filling in the antiferromagnet. 



■^°4{ v>.M|X)-Jh1 -'M^"'''''''' 

(AlO) 

remains small as long as |/i(x)| ^ Jh/U"^- When the 
latter inequality is violated, Eqs. (A9-A10) become in- 
valid. Using Eqs. (31), (32), and (A3), we write for the 
thermodynamic potential in the case of small xq , 



Via = -Ux^xi -dJ+—^n^~ 



(G) 




(All) 



The sum on the r. h. s. [which can be evaluated using 
the small-fc expansion and Eq. (A4)] is found to be of the 
order of xq( and can be omitted. Eqs. (A2) and (AlO) 
then yield the final expression. 
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It is easy to see that our neglecting the first term on 
the 1. h. s. of Eq. (A7) is appropriate as long as 
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At small absolute values of /x(x) < 0, this is satisfied as 

long as Xq ~ Jh ImI^^^^/^^ is small [see Eqs. (AlO) and 
(A4)]. On the other hand, when ij,{x) is close to fio {when 
Eq. (AlO) yields xg ~ ^x) - no\{J% + 4d2)/(L/d)}, 
Eq. (A13) takes form U > 7r(j|f + 16)^/2 j^ 2D and 
U[^l{x) - ^loY'^ > TT^{\Jl + 9)1/4 -^ 3o [gg]^ gjjjjjg ^jjg 
actual value of U /t for the CMR manganates is about 16, 
our results in the latter region can be viewed as an order 
of magnitude estimate only. 

As the value of x is increased towards half-filling, and 
the chemical potential reaches small (negative) values 

of /^(x) ^ —Jh/U"^, the value of xq starts increasing 
rapidly. The spin-up G-antiferromagnetic band is filled, 
XG = 1, for tJL{x) > dU/[2{JH + Uf] [cf. Eq. (32)]. In 
principle, the mean field equations can also be analyzed 
within this narrow area near jjl{x) = 0, assuming that 
Jh + Uxg is large. 
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FIG. 1. (a) Schematic representation of Eq.(9), with the 
two vertices corresponding to Tl'i and Ti.ii [see Eqs. (5) and 
(7)], respectively, (b) Second-order interaction correction to 
the magnon energy [cf. Eq. (18)]. (c) Second-order con- 
tribution to the temperature-dependent interaction F^^ be- 
tween two spin-up electrons, see Eq. (21). In all cases, solid 
and dashed lines correspond to electron and magnon Green's 
functions, respectively. Up- and down-arrows denote spin of 
the electrons. When evaluating these diagrams, one should 
ensure the proper antisymmetrisation of the the spin-up elec- 
tronic "legs" of each vertex, taking into account the appropri- 
ate momentum dependence. Momentum integration is greatly 
simplified in the large- Jg case considered here. 



FIG. 2. (a) Doping dependence of spin stiffness for a 
two-dimensional system with J = 0, Jn/t = 5 (solid 
and dashed lines, corresponding to U/t = 16 and U = 
cases, respectively). The dashed-dotted line corresponds to 
Jn/t — 14.6 and [7 = 0, and the dotted line - to Jh -^ oo. 
(b) The leading-order magnon energy in a 2D system with 
U/t — 16, J — 0, and x — 0.6. Solid, dashed and the up- 
per dotted lines are plotted using Eq. (19) and correspond 
to Jn/t — 5, Jn/t ~ 0.2, and Jh -^ oo, respectively. The 
dashed line is the result [43] for Jn/t = 5, U = 0, and the 
lower dotted line is the corresponding Heisenberg fit. 
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FIG. 3. Values of superexchange J corresponding to the in- 
stabilities of the ferromagnetic order at T = in a 2D system 
with Jh/I = 5, t/ = (a), and U/t = 16 (h). Dotted lines 
correspond to spin-wave instabilities, solid (dashed) lines - to 
the phase separation into G-type (A-type) antiferromagnetic 
and ferromagnetic phases, and the dashed-dotted line - to 
the phase separation into chain (see Fig. 4) and ferromag- 
netic phases. 



FIG. 4. Spin ordering in the chain phase. 
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